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Electronic excitations in metallic nanoparticles in the optical regime that have been of great 
importance in surface enhanced spectroscopy and emerging applications of molecular plasmonics, 
due to control and confinement of electromagnetic energy, may also be of potential to control the 
motion of nanoparticles and molecules. Here, we propose a concept for trapping polarizable particles 
and molecules using toroidal metallic nanoparticles. Specifically, gold nanorings are investigated for 
their scattering properties and field distribution to computationally show that the response of these 
optically resonant particles to incident photons permit the formation of a nanoscale trap when 
proper aspect ratio, photon wavelength and polarization are considered. However, interestingly 
the resonant plasmonic response of the nanoring is shown to be detrimental to the trap formation. 
The results are in good agreement with analytic calculations in the quasi-static limit within the 
first-order perturbation of the scalar electric potential. The possibility of extending the single 
nanoring trapping properties to two-dimensional arrays of nanorings is suggested by obtaining the 
field distribution of nanoring dimers and trimers. 


I. INTRODUCTION 

The development of Penning [I] and Paul traps 0 led 
to unprecedented measurements of fundamental proper¬ 
ties of charged particles by isolating individual electrons 
and ions Q. Singling out and having the ability to probe 
individual ions allows laser manipulation of their quan¬ 
tum states to form qubits Si , which is capitalized upon 
in the so called trapped ion quantum computer Q. Preci¬ 
sion measurements benefit from trapping or even slowing 
down the quantum systems, as reported by Veldhoven, 
et al., who carried out high resolution spectroscopy by 
slowing down a pulsed 15 ND 3 molecular beam in a Stark- 
decelerator Q. Deceleration can further help introduc¬ 
ing cold polar molecules into traps or storage rings Q. 
For a comprehensive account on the importance of polar 
molecules and the selection of those suited for decelera¬ 
tion and trapping experiments, see the recent review by 
van de Meerakker, et al. |;9(. Similarly, the development 
of optical trapping 0 of small particles, capitali zing 
on force gradients associated with radiation pressure [ll| 
exerted by tightly focused laser beams, led to the de¬ 
velopment of optical tweezers fl2j . Optical trapping in 
the tweezer configuration has produced significant vol¬ 
umes of data including those related to measurements of 
micro-rheological properties of aerosol particles [0 , op¬ 
tical binding of trapped particles [3] , and DNA and cell 
studies 0. 

The ability to manipulate the spatial distribution and 
time dependence of the electromagnetic fields interact¬ 
ing with charged particles or neutral molecules of fi¬ 
nite dipole moments is of prime importance in the fur- 


* msalhi@vols.utk.edu 
' passianan@ornl.gov 
t siopsis@tennessee.edu 



FIG. 1. The concept of a nanoring based toroidal trap. Cold 
polar molecules near the gray shaded surface approaching the 
central region may be trapped within a nanometer scale vol¬ 
ume due to the potential barrier established by the optical 
response of the torus to an incoming beam of photons with 
proper polarization, wavelength, and amplitude modulation. 


ther development of methods to contain and control sin¬ 
gle or ensembles of molecules jl6j. Local confinement 
and enhancement of electromagnetic fields, per recent 
advances in nano-optics and plasmonics 17j, will un¬ 
doubtedly prove important in the further development of 
systems that can exert electromagnetically induced spa¬ 
tial constraints on molecular and particulate populations. 
For applications in quantum information processing, it 
is important to have the ability to trap arrays of parti¬ 
cles jl8|. In particular, the development of arrays of elec¬ 
trodynamic nano-traps, as proposed here, may offer an 
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interesting perspective in storing and manipulating well 
separated but compact set of molecules. Advanced nano¬ 
fabrication 0 may be employed to generate a plethora 
of nanostructures [20( that can alter the electromagnetic 
fields to create specific force fields within which molecules 
may reside in slow motion. A variety of trap configura¬ 
tions and techniques have been proposed and designed, 
in particular for cold polar molecules pH - liUj . For on 
chip applications, electrostatic surface trapping of cold 
polar molecules have been proposed using various elec¬ 
trode configurations including a sing le charged circular 
wire [35|, parallel charged wires (36), and two charged 
rings above a grounded conductor plate embedded in an 
insulator f37j . 

In this article, we propose and investigate the con¬ 
cept of sub-wavelength electromagnetic traps in the form 
of metallic toroidal structures (micro- and nano-rings), 
shown in Fig. [l] Throughout this work, the term trap 
is defined to describe the electromagnetic environment 
in which a particle assumes a dynamic state with an 
eigenenergy spectrum insufficient to overcome a poten¬ 
tial barrier. Similarly, nano-trap refers to a physical ar¬ 
rangement in which appropriate particles experience a 
potential barrier with a spatial extension in the order 
of nanometers. We propose to achieve a nano-trap with 
toroidal geometry using noble metal nano-rings, which 
upon interaction with long wavelength linearly polarized 
photons produce a three-dimensional trapping potential 
barrier. The barrier shape is shown to be further con¬ 
trolled by use of circularly polarized photons. We envi¬ 
sion the toroidal nano-trap can be used to capture polar¬ 
izable dielectric nanoparticles, similar to optical tweez¬ 
ers 0, as well as cold polar molecules, assuming that 
the particles and molecules have been through a cooling 
process to reach a temperature of a few mK using an 
appropriate deceleration mechanism, for example by be¬ 
ing transmitted through a Stark decelcrator in the case 
of molecular beams as reported by Bethlem, et al. 39:). 
The studied structures are topologically distinct from the 
more commonly used geometries of a sphere, cylinder, or 
cube [H|: tori possess a distinct electromagnetic field dis¬ 
tribution that is here shown to allow for the formation 
of a trap for particles and polar molecules. As a com¬ 
parison, a quadrupole trap typically consists of a ring 
electrode with two end-caps. Here, in addition to the 
length scale difference, the proposed concept achieves an 
efficient trap in a configuration without end-caps. The 
designation toroidal nanotrap thus emphasizes the geom¬ 
etry of the structure employed to create a trapping region 
rather than the form of the trapping field. The actual 
trapping volume will be shown to take a rather spheroidal 
shape. Furthermore, comparing the proposed concept to 
quadrupole or hexapole electrode based traps (H, l3l| . 
it will be shown that, while the field lines in both cases 
may bear certain resemblance in the central region of the 
trap, the toroidal optically induced trap is solely due to 
the dielectric response of the metal nanoring. However, 
the advantages of optical excitation as opposed to apply¬ 


ing a voltage are manifold: a) optical trapping can be 
achieved off the chip or remotely without making phys¬ 
ical connections, b) optical fields can be provided with 
high modulation rates, polarization control, wavelength 
tuning, controllable power and intensity profiles, all of 
which can help achieving the desired field properties, c) 
optical fields can be delivered dynamically (as a rough 
comparison we may consider the method of laser scan¬ 
ning in digital optical disc data storage devices), d) mul¬ 
tiple lasers may be employed simultaneously, each with 
a desired polarization state, but with well defined phase 
relationship, and e) optical excitation is compatible with 
direct voltage application and thus in a proper configura¬ 
tion can be used in conjunction with the latter for other 
potential applications. 

The presentation is organized as follows. We begin by 
characterizing the optical response of a single torus in 
Section El In doing so, we will first treat the problem 
analytically in the quasi-static limit (where the field is 
generated by the scalar potential only) since the trapping 
properties are investigated in a spectral range where the 
wavelength of the force inducing photons will be larger 
than the dimensions of the nano-rings. Computational 
techniques are then invoked to obtain the time depen¬ 
dent scattering properties. By defining a trapping vol¬ 
ume in Section Ell we discuss the toroidal trapping en¬ 
ergy barrier for 15 ND, 3 , a candidate molecular system for 
deceleration studies |41| . Discussing the results in Sec¬ 
tion EO we also introduce multi-ring systems or arrays 
with the potential to control a lattice of polar molecules 
or particulates - an arrangement which may be useful in 
quantum computing and quantum information process¬ 
ing [lH and in sensing and imaging (42|. Concluding 
remarks are provided in Section [V] 


II. FIELD DISTRIBUTION AND SCATTERING 
PROPERTIES OF TOROIDAL PARTICLES 

The potential barrier experienced by a trapped parti¬ 
cle is assumed here to be the result of the interaction of 
the particle with an applied electric field modulated lo¬ 
cally by the torus. Since our objective is to describe the 
concept and potential use of the metal torus as a particle 
trap, we assume that electromagnetic properties of the 
trapped particle will not significantly affect the response 
of the trapping particle. This yields a decoupling of the 
problem such that the field of the trapping particle may 
be obtained while neglecting any modification induced by 
the radiative or polarizing effects of the trapped particle. 
The solid material domain making up the metallic torus 
is defined by a local dielectric function e (and conductiv¬ 
ity er). Within the spectral and sub-Kelvin temperature 
ranges considered, we will safely neglect the tempera¬ 
ture dependence of e, which otherwise has been shown 
to create thermoplasmonic effects in thin films [43] and 
nano-particles [44j |. To obtain the field modification en¬ 
gendered by the presence of the torus, we aim to char- 
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acterize the scattering properties of gold nano-rings, be¬ 
ginning with a single isolated nano-ring in vacuum. We 
therefore consider a torus, defined by fixing one of the co¬ 
ordinates u = wo in the toroidal system ( 77 , 77 , 0 ), which 
is related to Cartesian coordinates via: 

sy/u 2 — 1 sy/u 2 — 1 . s sin 77 

x = -cos (p,y = -sm 0,2 = - , 

u — cos 77 u — cos 77 u — cos 77 

where u £ [l,oo[, 77 ,0 £ [0,27r], and s is a scale factor. 
We first consider the response of the nano-particle in the 
quasi-static regime to an arbitrarily polarized uniform 
applied field E a = E a (sin 0, 0, cos0), where 0 is measured 
from the 2 axis. Introducing a perturbation parameter £, 
we expand the scalar potential \J/ a in the toroidal system 
as: 


= E a (x sin 9 z cos O') (1) 

OO 

= —‘X.ypi'K^ 1 sE a y/ w — £ cos 77 Gn(f , 9] £), ( 2 ) 

n—0 

where 

G n (r , 6\ 0 ) = sin 9 cos 0 cos 77,77 Q\- 1 («) 

+ cos 9 sin 77 . 77(1 — Sno)nQ° n _i(u), (3) 

with P™(u) and Q™(u) denoting the Legendre functions. 
Here, £ is introduced as a bookkeeping parameter for the 
perturbative calculations; at the end, we set £ = 1. The 
polarization 9 can be set by a polarization rotator and 
may be treated as a parameter so that 0 = 0 generates 
E a || 2 , while 9 = tt/2 generates E a _L 2 . The field E a is 
modified by the response of the torus to produce a total 
field, which we will study in two cases of polarization 
consistent with the toroidal symmetry. 


We now proceed to obtain the potential distribution in 
the case of E a || by superimposing tha and the solution of 
the Laplace equation. The potential of the trapping field 
E then takes the form: 


OO 

* =Ar E n sin nr] 

n= 1 


<2°iOo) 


i*_i(«) + Q°_i(«) 


(4) 

where Af(u) = — 2v / 27t l sE a y/u 1 and the coefficient a n 
can be obtained from the continuity of the potential and 
the w-component of D = eE across the toroidal boundary 
u = uq. At zeroth perturbative order (£ = 0), we obtain 
a„ = (1 - e)/(e - c n ), where 



L 2 J 



VuQ° n _ i(u) 



Similarly, in the case of perpendicular E± polarization, 
specifically E || x, the potential distribution of the trap¬ 


ping field is expanded as: 


4/ = A/"cos0 
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^ l(«)+Ql_±( U ) 


( 6 ) 


which at the zeroth perturbative order, produces al = 
(1 - e)/(e - cl) with 


SI ^g 

Vu p l i(u) 


Si l0 g 

VuQl 1 (u) 
n 2 



(7) 


In the quasi-static dipole approximation, the scattering 
cross section of the ring can be obtained by consider¬ 
ing only the dipolar contribution of the induced surface 
charge on the ring to the scattered field. Using the exact 
interior field and the exact exterior total field (as op¬ 
posed to the first order approximate field in Eq. ©), the 
dipolar charge density a c is found to be: 


( 7 , 3 ( 77 ,0) = (uq — SU 177) 3 / 2 ^ M n cos 77,77 cos 0, (8) 

n —0 


where 


M n = 


[A n -E a { 1 + M] 

7 r v 

Ql_i( u oY 


Ql- i(«o) - P n -l( u o) 


p Li( u o ) 


and the coefficients A n will be obtained from the conti¬ 
nuity of the displacement field D , yielding the following 
three-term difference equation: 


En+l QnE n T E n — 1 — An-t-1 2uqX n T Xn— 1, (9) 

where 


E n = 


o) - V n (uo) 


q n = 2u 0 + 


(e- l)Q^_i(«o) 

tQl-iM - p n{u 0 )' 


Xn. — 


Qn-1 ( u o) — r Pn(y-o) E a (l + 5 n o), 


( 10 ) 


which will be solved here using a Green function as de¬ 
scribed first by Love j4f|, and revisited by Kuyucak, et 
al. HU- Integrating the induced charge density over the 
volume of the ring using the surface charge, Eq. m yields 
the induced dipole moment, which is proportional to the 
polarizability a*. Thus, the dipolar scattering cross sec- 
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tion ads as given by: 


87TW 4 
a d s = ^r 


O-i 


( 11 ) 


can be computed. The absorbed power by the ring 
elevates its temperature T{f,t) which can be calcu¬ 
lated by solving the heat diffusion equation pC p T(f,t) — 
kAT(f,t) = s(r,t), where p 1 C p: and k are the density, 
specific heat, and thermal conductivity of the ring ma¬ 
terial, respectively. Any modulation information of the 
amplitude of the incident field will be embedded in s. 

The solutions above can be used to visualize the quasi¬ 
static toroidal fields and the dipolar scattering cross sec¬ 
tion within the zeroth perturbative order for the two 
polarization states of the incident photons, and to val¬ 
idate and further facilitate comparison to the computa¬ 
tional results for the more accurate time dependent re¬ 
sponse. Prior to obtaining the latter, we first employ 
the finite element (FE) method [47j to solve the Laplace 
equation in a three dimensional Cartesian computational 
domain, where an initial uniform static field is estab¬ 
lished by assigning equal and opposite potential values 
T| ±fc , k = xo,yo,zo to the appropriate domain bound¬ 
aries k corresponding to the two considered polarizations 
E a || and E a ±. We will thus solve V • (eV\I/) = 0 in 
the domain containing the torus by applying h ■ D = 0 
to all other boundaries. Following the quasi-static re¬ 
sults, to obtain the scattering properties of the single 
as well as several multi-ring configurations, we will solve 
V xH = d(eE)/dt+aE together with Vx£ = —dH/dt, 
employing the finite difference time domain (FDTD) 
method [Hj]. Here, again a Cartesian computational 
domain with absorbing boundary conditions (perfectly 
matched layers) is used to mimic the vacuum, and a lin¬ 
early polarized plane wave is propagated to interact with 
the torus. 

Having determined the response of the toroidal nanos¬ 
tructure, we will now describe how it may be suspended, 
and what the effect of such a suspension may be on 
the electric field of the nanotrap. For any experimen¬ 
tal attempts to trap molecules, the nano-traps need to 
be suspended without significantly distorting the electric 
field of the trapping region. As a preliminary suspension 
method, one may consider a membrane containing single 
or arrays of toroidal traps. Such an array of nano-traps 
may be fabricated using a low loss substrate. Employ¬ 
ing a focused ion beam lithography approach, one may 
attempt to carve channels in the membrane beneath the 
opening region of the torus, or fabricate the traps on 
other appropriate meshed substrates. For an optically 
transparent membrane, we anticipate negligible scatter¬ 
ing or modification of the local electric field. Studying 
the optical responses of metallic nanoparticles, similar 
configurations arise when using a photon scanning tun¬ 
neling microscope (PSTM) or a near-field optical micro¬ 
scope (NSOM) for imaging purposes, or for mapping the 
local field distribution. In such cases, when an optically 


conductive nanometer sized probe is brought in the prox¬ 
imity of the metal nanostructure (in our case, the torus) 
to register the field distribution, the question of what ef¬ 
fects the approaching probe or the substrate on which 
the nanoparticle resides may have on the nanoparticle 
surface modes has arisen in recent investigations. Such 
effects can be shown to be negligible for non- or low- 
absorbing substrates, such as fused silica, glass, and for 
parts of the spectra, silicon and germanium. Although 
the presented concept is based on an operation away from 
the plasmon band of the nanostructure and without a 
strong need for spectral tunability, in plasmon related 
studies, the (dielectric) substrate effect on the plasmon 
energy spectrum typically only amounts to a minor red- 
shift. Therefore, substrate effects can be safely assumed 
to be small, and more importantly, non-detrimental to 
the trap formation. More specifically, with reference to 
Eq. 3 in {4f|, a similar study can be performed to quan¬ 
tify the negligible effect of a low-loss substrate. 


III. ELECTROMAGNETIC TRAPPING WITH 
NANO-RINGS 

Our objective here is to investigate the potential of 
metal nano-rings in applications related to electromag- 
netically trapped species possessing a permanent or in¬ 
duced electric dipole moment ft. We wish to show that 
microscopic objects carrying a moment will be unable to 
escape from a trapping volume in the vicinity of the cen¬ 
tral region of the torus. We assume that the particles 
to be trapped are in thermal equilibrium such that their 
positions and velocities follow an appropriate distribu¬ 
tion. Nanoparticles are readily aerosolized at room tem¬ 
perature, however, in this study, the assumption is that 
they, too, have been cooled down and introduced in the 
toroidal region. Our primary focus here is on trapping 
molecules. The field E distribution of the ring is thus 
to present a potential energy barrier to the molecules. 
The engendered field by the metal nano-ring perturbs a 
molecule in a rotational state with an interaction energy 
of the form He = —ft ■ E modifying the rotational en¬ 
ergy levels of the molecules (Stark perturbation). For 
certain polar molecules used in Stark deceleration stud¬ 
ies, the Stark energies of the rotational states have been 
calculated as a function of the applied field strength Q . 
Specifically, for our objective, the assumptions regarding 
the quantum states of the molecules are similar to those 
leading to the Stark energy expression for ammonia iso- 
topologue 15 ND 3 (see Eq. (11) in ref. Q): 


Wstark(u,»7>0) 



( 12 ) 


where T can be taken from Eq. [6] or be provided com¬ 
putationally from FE and FDTD, W; nv denotes the en¬ 
ergy splitting associated with the atomic inversion within 
the polar molecule in the state \J,K,M), and a = 
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FIG. 2. Polarization dependence of the field at the center of 
the torus. For a gold nano-ring of appropriate aspect ratio 
7 , when the polarization of the incident held is perpendicu¬ 
lar to the symmetry axis of the torus, the electric held can 
be made negligibly small. The analytical solution using hrst 
order perturbation, while predicting a lower value for the as¬ 
pect ratio for this minimal held, overall it agrees well with the 
computational solution employing the hnite elements method. 


J( J +1 )/yKM is related to the matrix element diagonal 
in K and M of the Stark perturbation such that the first 
order Stark energy is —E/a = (J, K, M\ He \J, K, M). 
Similarly, second-order Stark energy correction and high 
field strength Stark effect including higher-order pertur¬ 
bation terms may be considered as necessary. Using 
Eq. [12] and the held profile of a nanoring with given ge¬ 
ometric characteristics, the energy barrier along the x, y 
and z directions may be studied. Furthermore, for the 
sake of future studies, and without further elaboration 
here, denoting the mass of the ND 3 molecule with m, 
forming the Lagrangian and obtaining the equations of 
motion: 


(■u 2 


ms 2 u 

l)(u — COST ]) 2 
ms 2 ?) 

(it — cos r \) 2 
ms 2 (u 2 — l)</> 

{u — COST ]) 2 



cWstark 

= 0 

(13) 

+ 

du 


cWstark 

= 0 

(14) 

+ 

drj 


<9Hstark 

= 0 , 

(15) 

+ 

8(j) 


the trapped molecular motion may be studied within the 
trap volume or molecular deceleration may be investi¬ 
gated for molecules moving through the trapping region. 


IV. RESULTS AND DISCUSSION 


We begin our study of the trapping properties of metal 
tori by characterizing the optical response of single gold 
ring tori in the spectral range from visible to the short 
wavelength infrared. We will then extend the results 
by considering more elaborate configurations of dimers 



FIG. 3. Comparison of the quasi static and time dependent 
electric field across the torus. For an aspect ratio 7 = 3.5, 
the analytical results as well as the computational results in 
both the electrostatic (FE) and transient (FDTD) limits con¬ 
sistently produce a field profile suggestive of a trap forma¬ 
tion. The held lines computed using the FDTD method are 
noted to display smooth curvature, whereas both the analytic 
and quasi-static computational approaches produce minor ar¬ 
tifacts around the low hied regions of the trap. 


and trimers as building blocks for arrays of nanorings for 
both trapping and deceleration of cold polar molecules. 
To study the geometric dependence of the trapping po¬ 
tential, we introduce 7 = R/r = uq, the aspect ra¬ 
tio, where R = auo/y/u 2 —l is the major radius and 
r = a/ \Jv/q - 1 is the cross sectional radius. The material 
making up the torus is taken to be gold, although other 
suitable materials may be considered. To determine the 
field distribution and energy dissipation by the proposed 
structures, we will resort to simulations, based on the an¬ 
alytical results, and computation employing the FE and 
the FDTD methods. For the latter, the Drude model 
is numerically adjusted to fit the values of the complex 
dielectric function of gold, e, reported by Johnson and 
Christy [50|. The fitting parameters ujj = 2155.6 THz, 
7 e> = 18.36 THz and too = 9-06 reproduce the values in 
the spectral range 0.4 — 2.2 y,m, which will be considered 
for the FDTD simulation of the electric field and absorp¬ 
tion spectra of the various ring configurations studied. 

The results of both the FE computation and the per¬ 
turbative analytic calculation for the held value at the 
center of the ring and its dependence on the aspect ratio 
are plotted in Fig. [2] The results, being in good agree¬ 
ment, suggest a special value of the aspect ratio (7 ~ 3) 
in the case of perpendicular polarization of the incident 
beam ( E a J_ z) for which the held is almost extinct. Fur¬ 
ther analysis sets this values closer to 7 « 3.5, for which 
we then compare the results of the FE, FDTD and per- 
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FIG. 4. Computed scattering (upper panel) and absorption 
(lower panel) cross sections for single gold nanorings of dif¬ 
ferent aspect ratios 7 . For a fixed minor radius r, the larger 
aspect ratios exhibiting the larger cross sections correspond 
to rings with larger main radii R. 

turbative analytic calculations by plotting the field along 
the polarization direction for each method in Fig. [3] Ev¬ 
idently, for the geometric configuration at the given as¬ 
pect ratio, the field distribution suggests the possibility 
of forming a trap. 

In the time domain, the interaction of either a contin¬ 
uous wave (CW) or an amplitude modulated (AM) laser 
beam with the torus results in a total absorbed power P t . 
Assuming that a linearly polarized field E(r,t ) of wave¬ 
length A and propagation vector k can be provided by a 
an appropriate laser, we specify the incident radiation to 
be in the form of a beam of a circular cross sectional area 
A, pulse duration At, the repetition rate T) = 1///, en¬ 
ergy per pulse e p and an average power P „. For a single 
gold nanoring with an aspect ratio 7 and volume V, de¬ 
noting the incident intensity the absorption and scat- 



FIG. 5. Field profiles along x, y, and 2 directions across the 
torus for various aspect ratios for an incident field Ei polarized 
perpendicular to the symmetry axis 2 of the torus. The insets 
visualize the field distributions in the corresponding planes 
for a ring with 7 = 3.5 interacting with an incident field 
of a wavelength A = 2 ym propagating along 2 . The field 
distributions were obtained from FDTD computation. 
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FIG. 6 . (a) Stark energy profiles along x, y, and z directions 
through the center of the torus corresponding to the energy of 
ND3 molecules interacting with the total toroidal field, that 
is, the incident field Ei polarized in the plane of the torus and 
the scattered field. The field distribution was obtained from 
FDTD computation, (b) Stark energy profiles through the 
center of an vertical dimer with constituent rings separated a 
distance d = 110 nm and having an aspect ratio 7 = 3.5. 


tering cross sections are defined as a a bs = J v pdV/Ii and 
&scat = f s s 'scat ■ dS / R, respectively, where p is the ab¬ 
sorbed power density in the ring and S sca t is the scattered 
Poynting vector over its surface. Neglecting photolumi¬ 
nescence, the extinction cross section is a e = a a bs + Cscat 
and the total absorbed power by the particle is given by 
Pt = tJabsIi ■ Obtaining the scattering and absorption 
cross sections is relevant to our discussion, as it com¬ 
pletes our understanding of the optical response of the 
nanoring to an incident wave. The response is critical to 


the feasibility of the nano-trap. Scattering and absorp¬ 
tion cross-sections are plotted for different aspect ratios 
and wavelength range A £ [0.4, 2.2] p ,m in Fig. 2] The 
electromagnetic trap that we are discussing is illustrated 
by plotting the field profile in the vicinity of the nanoring 
along the x, y and z directions for different aspect ratios 
in Fig. [5] 

To describe the trapping of cold polar molecules, we 
will consider ND3, where the molecular beam is in a sin¬ 
gle weak field seeking state | J. K, M) = 11,1, — 1). The 
ND3 system is known to have an electric dipole moment 
of p = 1.48 D and a frequency W- lnv = 1.43 GHz at which 
the position of the nitrogen atom inverts [ 9 }). To min¬ 
imize escape from the trapping potential, the molecular 
population can be prepared to have appropriate mean- 
free-path so as to render molecule-molecule encounters 
less likely. For the trapping volumes considered here, 
this can be readily achieved already in the low vacuum 
regime. Assuming that when the applied electric field is 
switched off, that is, when E = 0, the energy of a free 
cold molecule of mass m and temperature T can be sam¬ 
pled from a Maxwellian distribution, it will have a root 
mean squared speed of v = (SksT/m) 1 ^ 2 m/s, where fee 
is Boltzmann’s constant. For T = 25 mK [21j a velocity 
of v = 6 ms -1 implies that a cold molecule may depart 
from the trapping region of a ring with 7 = 3.5 and 
R = 87.5 nm. When E ^ 0, a threshold trapping energy 
W s trk in Eq. (fT^l) corresponds to a minimum field mag¬ 
nitude of E th = atVF^k/W^k+Wnv)] 1 / 2 . From decel¬ 
eration studies by Bethlem, et al. [2ll| . Wg^ ark > 25 mK, 
yielding a minimum field amplitude of E th > 0.27 V//im. 
Such a value for the local electric field, which can be es¬ 
timated to be of the same order of magnitude as the 
incident field due to the spectral position of the plasmon 
band being far from the laser wavelength, can be made 
available from an appropriate continuous wave (CW) 
laser. We also note that high resolution infrared ab¬ 
sorption spectra of 15 ND 3 has been reported by Cane, 
et al. f5l| and thus the incident laser wavelength can 
be tuned accordingly so as to minimize excitation of the 
molecular bands. 

Furthermore, for a given field amplitude E th , the total 
absorbed power P t = i cneo<j a b s \E th \ 2 will be below the 
melting point of the nanostructure. Since for the typical 
rings considered here the volume is V ~ 1.1 x 10“ 15 cm 3 , 
we assume that the deposited energy is distributed uni¬ 
formly, i.e., s(f,t ) = s(t) = PR is the source term ob¬ 
tained from the total absorbed power by the nanostruc¬ 
ture. Employing FE method, we here computationally 
verified the temperature uniformity by solving for the 
time dependent thermal response of a gold nanoring with 
a heat capacity C p = 129 J K _1 kg _1 , a thermal conduc¬ 
tivity k = 317 Wm _1 K _1 , and a density p = 19300 kg 
m . The results also predict negligible temperature in¬ 
crease due to the absorbed power of a silicon substrate- 
supported gold ring in vacuum but otherwise subject to 
mK near substrate boundary temperature. 

Using Eq. m, the stark energy and the field profile of 
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FIG. 7. Visualizing trapping boundaries for the 15 ND 3 system 
under circular polarization excitation. The field isosurface 
rendering a single ring trapping volume (left) of 2.7 x 10 s nm 3 
and a ring dimer trapping volume (right) of approximately 
5.5 x 10 s nm 3 . At these boundaries, molecules in a single 
weak field seeking state \J,K,M) = 11,1, — 1) would possess 
a Stark energy of Wstark = 10 mK. 


nano-rings with an aspect ratio 7 = 3.5, we plot the cor¬ 
responding energy barrier along the x , y and z directions 
for two systems, with the choice of the aspect ratio based 
on the depth of the potential energy barrier. The first 
system is a single gold nano-ring, whereas the second sys¬ 
tem is made of two axially aligned gold nano-rings with 
an optimized separation distance d = 110 nm for higher 
energy barrier as illustrated in Fig. [ 6 ] The dimer ring 
trapping system has a higher energy barrier when com¬ 
pared to the single ring system. The latter is a realistic 
trap for molecules with a translational energy along the z 
direction of the order of few mK. The dimer trap frequen¬ 
cies reaching the GHz range may contain several tens of 
bound levels (about 75 at about 1 GHz [ 52 }). To compare 
the actual trapping volume in both cases, we note that 
both the position and shape of the trapping boundary 
can be affected by the specific lTg^ ark and the incident 
photon polarization, as shown in Fig. 0 for the case of 
circular polarization. 

The proposed concept can be extended to include an¬ 
other important case where the toroidal nano-trap may 
be considered as a potential building block of a quantum 
computer, which requires the manipulation of qubits, as 
discussed, e.g., by DiVincenzo [53|. The most elusive re¬ 
quirement for building a quantum computer is to have 
a set of individually addressable qubits. This translates 
into the formation of a 2D array of nano-rings. The array 
is characterized by a minimum separation distance that 
would preserve the electromagnetic trap profile, found in 
our example simulation to be d m i n ~ 0.5 /im between 
individual centers of the rings. Thus, to verify that the 
individual rings are not coupled in a way that would alter 
the trap profile, we obtain the absorption and scattering 
cross section for an array of rings with a minimum separa¬ 
tion d m i n . Both plots are shown in Fig. [8] One can easily 



X [nm] 


CM 

E 


o 

c/) 


o 

k_ 

o 


3.5 


2.2 



FIG. 8 . (a) Computed scattering and absorption cross sec¬ 
tions for single planar gold nanoring trimer with each ring 
having 7 . With the rings sufficiently apart, the cross sections 
approach those of a single ring, (b) Comparison of the ab¬ 
sorption and scattering cross sections for a vertical dimer with 
the same parameters as those in Fig. [5] 


see that they are almost identical to the absorption and 
scattering plots obtained in the case of individual rings. 

In closing, we note that an experimental exploration 
of the proposed concept may begin in a configuration 
of scanning probe microscopy. For example, the hybrid 
photonic-nanomechanical force microscopy (HPFM) or 
the near-field scanning optical microscope (NSOM), both 
with spatial resolutions below the diffraction limit and 
allowing high levels of control and manipulation of light 
at the nanoscale, may be employed to first characterize 
the electromagnetic response and morphology of single 
nanorings. Similar modalities permit spectroscopic in- 
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terrogation by both delivering the excitation photons lo¬ 
cally to the trapped molecules as well as collecting the 
emitted photons by the molecules. 

V. CONCLUSION 

The theoretically obtained scattering properties and 
field distributions of toroidal gold nano-particles strongly 
suggest that trapping of cold polar molecules and dielec¬ 
tric beads may be feasible optically in the proposed ge¬ 
ometry. For a wavelength of 2 /im, the results for a spe¬ 
cific value of the toroidal aspect ratio indicate the pos¬ 
sibility of a strong and well defined trap for cold polar 
molecules when the external field polarization is parallel 
to the axis of the ring. Interestingly, the results show 
that a trap forms away from plasmonic resonances of 
the nano-ring. The frequency dependence of the nano¬ 
ring response shows an enhancement of the electric field 
near the plasmon resonance band, as may be expected. 
However, at resonance the trap becomes shallow. Fur¬ 
thermore, from the single torus results, we suggest more 
elaborate configurations in the form of two-dimensional 
arrays of nano-rings toward the potential implementation 
of multi-particle nano-traps. Such compact and minia¬ 


turized structures may be of importance in the design of 
future hardware platforms for quantum computation ap¬ 
plications. Considering the large forces on the molecules 
achieved in the nanotrap with trap frequencies up to a 
GHz, it is conceivable to resolve these levels spectroscop¬ 
ically. The utility of the toroidal trapping concept may 
be extended to confining dielectric particulates via force 
gradients similar to dielectrophoresis [54j |. by creating a 
toroidal potential energy barrier. Assuming a particle 
with no internal degrees of freedom but with a dielec¬ 
tric function e p , the force leading to stable small ampli¬ 
tude oscillations of the particle may be obtained from 
Maxwell’s stress tensor T by integration over the parti¬ 
cle surface f(t) = ff T • ndA. Thus, particles such as 
polystyrene nano-spheres may be trapped via a similar 
approach. 
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